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Abstract 

Three body heavy meson decays allows one to extract important quantities for 
flavour physic, however the phenomenological analysis has some open problems to 
be solved in order to make possible to extract physical parameters with precision. 
After many years and some experience accumulated, we have learned new features 
and we expect to learn much more from data coming mainly from B three body 
decays. In this paper we address these aspects and we present a proposal to modify 
the Dalitz fit amplitude method in order to allow further studies. 



1 Introduction 

Three body heavy meson decay analysis have been treated as a coherent sum of indi- 
vidual two body resonances amplitudes plus a non-resonant component, each amplitude 
with a strong final state interaction phase (FSI). This simple configuration based in the 
dominance of the quasi two body component, was proposed in the early eighties by Mark 
collaboration [1]. Important changes were made to improve this approach during these 
already three decades, however remains the basic structure of the original model. This 
analysis method have been working well to represent low statistics charm three body 
decay but have been presenting some problems to represent the equivalent three body 
beauty decays. 

The basic concept used by Mark collaboration in its initial three body analysis was 
inspired in the amplitude analysis used for vr - nucleon and K - nucleon interaction giving 
three or four final state [21 S H El El E] B Nevertheless, there is an important difference 
between hadron hadron interaction and heavy meson decay: while the last one use a 
full coherence sum, in hadronic collisions this quantum effect had been used in a more 
restrictively form, assuming partially the coherence between different two body resonances 
ampUtudes present in the same phase space region. 

Recently, for some charm meson three body decay analysis [9l[Tn], a significative change 
was partially introduced by FOCUS collaboration. To take in account interferences among 
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resonances with the same quantum number, FOCUS used the two body unitary K-matrix 
approach to represent the s-wave component. 

The K-matrix in Dahts analysis and the usual approach of coherent sum of resonances, 
reflected two different treatment for particle interferences well characterized by Azimov in 
a recent paper [H]. According to this author, interference between particles with the same 
quantum number would be classified as "direct interference", otherwise it would said to be 
by "re-scattering interference". Among other differences between they two, the locality 
is one of more important characteristic pointed out in this paper. Direct interferences 
have a large phase space range, while re-scattering interference between two resonances 
is limited at the region where resonances have intercepting phase space. 

However, if for one side the K-matrix approach tried to treat correctly interferences 
between resonances belonging to the same partial wave, on the other side neglecting 
possible re-scattering interferences between particles with different quantum numbers, 
that is a singularity triangle type diagram [12j. This approximation seems to work for 
Df — 7- 7r^7r+7r+ decay [13j, but not for K^-k^-k^. In fact the Kt^ phase motion 

of the last decay, extracted with Partial Wave Amplitude (PWA) analysis first provided 
by E791 collaboration and then by CLEO and FOCUS [HI [TSl [IS], shown a distribution 
substantially different from the one observed in scalar partial wave Kir elastic scattering 
by LASS collaboration |T7] and also different from the recently PWA analysis performed 
for the semi-leptonic K'n^e'^Ue decays [18] %. However the same analysis for 

the 7r^7r+7r+ decays seems to have similar phase motion behaviour than the one 

obtained by 7r^7r+ elastic scattering [13] . 

Since both kinds of interferences, with different behaviour, can be present in three body 
heavy meson decay, a complete amplitude analysis has to take in account the existence 
of both effects and also provide from the fit parameters, the relative contribution of each 
effect for each intermediary state. 

2 Isobaric model and experimental results. 

With the increase of statistics, more resonances with different spins started to contribute 
significantly in three body D decays, Argus collaboration introduced a generalization of 
the isobar model [20]: each possible resonance amplitude is represented by a Breit-Wigner 
multiplied by angular distribution function associated with the spin of the resonance. The 
Non- Resonant is parametrized as a constant. The various contributions are combined in 
a coherent sum with complex coefficients that are extracted from maximum likelihood fits 
to the data. The absolute value of the coefficients are related to the relative fraction of 
each contribution and the phases takes into account final state interaction (FSI) between 
the resonance and the bachelor pion. This phase is considered constant because it depends 
only on the total energy of the system, i.e. the heavy meson mass. 

With this simple model plus the E791 method [2T] to determine the mass and the 

^ Both, the semi-letonic and the Ktt elastic scattering scalar phase motion agree each other in a clear 
verification of the Watson theorem [TH] 
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width of some not yet well defined resonances, most of three body charm meson decay 
were successfully described 0. On the other hand, B charmless three body decay scenario 
is not successfully describe by the use of the Isobar model. The most important problem 
is associated with the non-resonant (NR) contribution. Differently from charm decays, 
the NR contribution has playing a significant role in B meson studies. At a very early 
stage of the B factory analysis programme the 'constant amplitude' description of the NR 
term was abandoned and replaced by alternative approaches. When studying the decay 
B~^ K^7^^7r~ the Belle collaboration proposed an empirical exponential function to 
represent the NR contribution [2l|. BABAR instead, used for the same decay mode |25j a 
phase and magnitude parametrization obtained by the LASS experiment [17] for S-wave 
Kti elastic scattering. 

The empirical function proposed by Belle |23] has two terms, one for each Dalitz plot 
variable (sis = M'^{K~^Tr~) and S23 = M'^{n~^n~)), and is given by: 

Anr{K+n+n-) = <"e-""^^e<'' + a^^e-'^^^e'^^'' , (1) 

with five free parameters: one magnitude and one phase for each of these terms and 
the slope of the exponential function a. The fit results exhibit three noteworthy and 
correlated characteristics: i) a large NR fraction of around 34%, ii) a phase difference 
between 6^^ and 62'' around 180°, iii) a sum of fractions from all components, resonant 
plus NR of around 150%. In an amplitude analysis the sum of fractions can of course be 
bigger or less than 100% and this is not inconsistent with the unitary concept 0. However, 
such an outcome indicates the existence of destructive or constructive interference among 
two or more amplitudes. In practice, a sum of fit fractions ^ 100%, that is a large 
destructive interference, is often indicative of underlying problems in the fit model. In 
the study under consideration it may be that the apparent large sum of fractions is in 
fact an artifice, and that in the coherent sum the significant NR contribution is largely 
"eliminated" through the destructive interference arising from the ~ 180° phase difference 
between the two NR amplitudes. 

Subsequently, several amplitude analyses, including Belle's study of the mode B^ 
K'^TT+n- [2B] and the BABAR analysis of 5° KsK+K' [27], have used this exponential 
amplitude. The fit results for share the same three characteristics discussed above for 
B+ ir+vr+vr-. 

This approach has been also used by both BABAR and Belle for three-body decays 
involving two identical particles, namely B^ K^K^K~ [2U |28] and B^ — ^ tt+tt+tt^ 
|29] 1^. Again, in both analyses a large contribution is found for the NR component, and 
the overall sums of fractions are significantly above 100%. 

There is at least one important exception, the K^tt~tt~^. Because CP violation studies this 

decay have been analysed by the B-factories with a large samples and also many problems [22l [23] 

^ The concept of fraction in a Dalitz plot analysis is directly related to the branching ratio of a 
particular amplitude. It can be classically associated to that fraction of the total number of events in an 
incoherent sum of amplitudes. 

^ In these cases the expression Eq. ^ contracts to a single amplitude because of the symmetry under 
the exchange of identical particles. 
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Another important problem for almost all the charmless three-body B decays is that 
it has been necessary to include some scalar contributions to get a good fit to data. 
The mass and width values obtained in some analysis involving 7r"^7r~ in final state do 
not match among these two experiments [25] and some times do not match between 
different channels of the same experiment [29]. The same situation occurs with high mass 
scalar contribution involving the final state K^K~ [28]. 

The non-resonant component which, in general, spread within all the populated phase 
space, can mimic other dynamical components |30] through the interferences with the 
resonances present in the same phase space. High mass scalar resonances have also similar 
capability, due this angular distribution and the large width, these contribution populates 
homogeneously large portion of the di-hadrons phase space, without a clear signature like 
vector or tensor contributions. This shadowing phenomenon was observed in charm meson 
three body decays at the E791 experiment, where the overestimated contribution of the 
non-resonant amplitude use to be able to mimic the real contributions coming from the 
scalar mesons cr and fi:[2Tl [3T] . 

3 Interference among resonances with the same quan- 
tum numbers. 

In the isobar model, regarding interferences, all contributions are treated in the same way, 
independent of its origin. On the other hand K-matrix approach proposed by FOCUS 
[9l [10] use a especial treatment only to represents the S-wave scalar amplitude contri- 
bution unitary amplitude extract from two body elastic scattering. This approach take 
into account the peculiarity of long range phase space interference between two body am- 
plitudes with the same quantum numbers. The high spin resonances was treated as the 
isobar model. 

There are alternative ways to write two or more resonances with same quantum num- 
bers in a unitary fashion, Tornqvist [32], per example, proposed to use a couple Breit- 
Wigner. There is a more recent representation that can be very useful for three body 
amplitude analysis, it was proposed by Svec [33] based in a Hu method [51]. As we can 
see below, this approach can be useful, because it makes explicit the unitary constraint in 
the isobar model, while keeping its usual structure. The proposal consist to represent a 
multi-resonance contribution through a unitary product of isolated Breit-Wigner. Since 
one resonance is represented by an unitary S-matrix Breit-Wigner function the product of 
Breit-Wigner will be unitary also. This idea is suitable to Breit-Wigner or other amplitude 
to represent a resonance. 

With this assumption one can arrive to a simples T equation for two resonances [33] : 




where BWi is the usual Breit-Wigner function. This equation can be re-written as co- 
herent sum of Breit-Wigners with specific complex coefficient, depending of the resonances 




T 

■J- r 



res 



BWi{s) + BW2{s) + 2iBWi{s)BW2{s) 



(2) 
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parameters as we can see below: 



Tres = Cr>{s)BW,{s) + C^^'>{s)BW2{s) (3) 

the coefficients are given by Ci^\s) = 1 — 2i / {zi — Z2) and C'2\s) = 1 + 2i/{zi — Z2), 
where Zi = mf — imiTils) with a smooth dependence in the energy variable. So to 
guarantee unitarity, the phase differences between two resonances with the same quantum 
number should be write as: 



cP(.)-cj''(.)= ^'t'"?!:it"' w 

where m is the mass of the resonance, r(s) the width, Am^ and AmT{s) are respec- 
tively the mass and the mT{s) differences between the two resonance. Since the Am? is 
in general bigger than the AmT{s), and the only dependence in s variable comes from the 
width, one can expect a little changes in the phase difference. In some sense. Equation 
2 justify the use of coherent sum of Breit Wigner with constant parameters, at least for 
charm three body decay taking in account the reduct phase space of this decay. With the 
parameters extracted from data, one can check whether or not these values are compatible 
with expected by a two body unitary behaviour of resonances in the same partial wave. 

The suitability of the this approach to represent data involving amplitudes with the 
same quantum numbers, could be observed in two occasions: the scalar amplitude and 
the phase motion obtained with the partial analysis analysis (PWA), performed for the 
— )■ K'TT^TT^ [m [16] and for — )■ 7r~7r^7r+ [13] decays. Both decays were 
well represented with a coherent sum of two scalar resonances contributions with a phase 
difference, as it was plotted in references [131 [IH 113 HB]- However the K-matrix was able 
to represent well only the D"^ — ^ 7r^7r+7r+. The scalar K^-k^ amplitude and the phase 
behaviour of the — > K^-k^-k^ decay, obtained with PWA differed substantially from 
the i^~7r+ elastic scattering experimental results [T7]. This observed difference should be 
due a strong re-scattering contribution in this three body final state [35] , that can not be 
represented in the K-matrix formalism [12] . 



4 Interference among resonances with different quan- 
tum numbers. 

We can divide the interference between resonances with different quantum numbers in 
Dalitz plot in two types: parallel, i.e. defined at the same Dalitz variable Sij and that 
can interfere in a large portion of the phase space, like the vector resonance p(770) and 
the scalar one (t(500) or with the /o(980) in a decay involving two charged mesons tt in 
final state. The parallel interference have been described in several models through local 
interferences for all phase space, but due to the angular distribution, the phase space 
integral should close to zero to parallel resonances with different spin. 
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The second type we can call crossing, i.e. interference between resonances that cross 
each other in the phase space with different Sij Dalitz plot variables, like the interference 
between the p(770) resonance with the i^'*(890) in decay in K~^n^7r~. The crossing 
resonances has the interference in a limited phase space region and the integral in general 
is different from zero. With a coherent sum, this term make a real contribution to the to- 
tal amplitude square. Due to the formal complexity to work with this contribution, some 
Isobar approach to describe hadron hadron interaction, preferred to make the approxima- 
tion that interferences between crossing resonances zero. This approximation was used in 
SLAC manual for three body partial analysis in the eighties [7]. To performed the analysis 
for the final state K~7c~^tt~p, they assume the integral of the crossing term between the 
p and K*{890) is zero, beside they can observe a significant amount of interference. 

Since the crossing interference can not be estimated by first principles, a correct in- 
clusion of the crossing contribution in amplitude analysis is a big challenger. The best 
known formal approach to use for study re-scattering effects is the Faddeev equation 
[TTj |36] . Very preliminary studies using this method, have presented nice results with 
a natural explanation for the observed phenomena to a global phase difference between 
the P and S wave from Kit partial wave analysis (PWA) of the — )■ K'tt^tt^ decay 
to the Kn elastic scattering [37J. In fact other than the difference pointed out before 
between these two experimental results, there is a global phase difference bigger than 7r/2 
degrees[Tll |T5l [16]. This difference between three body decay and elastic scattering was 
also observed in h^h^ti^ partial wave analysis [T3] . 

Experimental studies in the sixties and seventies presented a practical and empiri- 
cal approach for the re-scattering problem. They treated this effect by a coherent sum 
of Breit Wigner also to the crossing resonances. In K~^p K^n^p studies at CERN 
experiment [3], they introduced a sum of individual amplitudes with a constant phase 
difference, between the interfering resonances due a re-scattering phenomena. This deter- 
mine the degree of constructive or destructive interference among them. Each amplitude 
is parametrized by a Breit- Wigner multiply by a real density term independent of the en- 
ergies, that can be extracted from data. However there is an important difference between 
the original amplitude analysis and the one we use nowadays, the coherence sum was not 
generalized and fixed for all resonances, like we have been doing in the three body heavy 
mesons decays. 

In fact in their studies, the coherence is include in the fit function only among reso- 
nances that they observe a priory in data, a signature of interference between them. For 
the K~^p K^TT^p study, they use a coherent sum only for the K*{1400)p and a baryonic 
resonance amplitude i^'°A^*(1236), but not between the dominant contributions K*{890)p 
and K^N* {1236). This approach was generalized by the Berkeley group that studied the 
interaction tt~p tt^tt^p [1]. In fact they introduced a coherence factor to take in ac- 
count the observed re-scattering for all resonances present in the spectrum. The observed 
results between the p~'s amplitudes and the different A^* resonances variety from 0.6 to 
0.8. 

Other physics process involving particle interference among different quantum numbers 
used also a free parameter to manager with the non obvious source of coherence. In 
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order to observe the p and u isospin violation interaction proposed by Glashow [55] . 
experimental studies were performed with low transverse momentum events from TT~p — )■ 
TT'TT+n interaction [39]. The amplitude square used to fit the 7r"'"7r~ spectrum, included 
one Breit-Wigner square for the p decays and another one for the G-parity violating 
channel uj tt^tt^, plus an interference term. The last one was wrote to represent the 
clear interference observed in the p mass region. Besides the usual product of two Breit- 
Wigners and a relative strong phase, was also included a real multiplicative coherent 
factor that could change from to 1. This factor was also include in other studies p and 
oj interactions studies with different reactions [^ HT| H2]. 

More recently was proposed also a coherence parameter to take in account possible 
contamination in the dominant pseudoscarlar exchange process, from the axial meson 
component in the ir'p vr'^vr^n interaction [33] . Also Cristal Barrel collaboration [H| H5] 
used a coherence parameter to take in account the many different partial waves in the 
initial and final pp interaction at rest, that part interfere and part does not. For the 
pp — )■ TT^rfrj [l5] they fit seven intermediary states, with seven complex number for each 
one resonance contribution plus twelve coherence factor. Some of them were set to zero 
because there was no overlapping region or apriori tests did not present any significant 
contribution. 

Heavy meson three body decays has a simpler configuration since that it involves an 
spin initial state and three spin zero final state. So the use of a coherence factor is 
less crucial than used in the above examples, were some times it is associated with and 
average of the initial or final helicity state for interaction involving baryons. However some 
problems remain in heavy meson decays, like the re-scattering, the existence of different 
isospin amplitudes in a same final state. 

Also a technical problem can coming out in amplitude analysis. In fact, experimental 
resolution have not been considered fully in their analysis. In general it is not important 
to describe the mass spectrum of the usual resonances, as it was pointed out for the 
narrow /o(980) resonance analysed in the — )■ h'^tt~'k~^ [IB]. However it might be 
different for describe the crossing region between two narrow resonances, the interference 
can created a strong variation in a short little phase space region and the experimental 
resolution can spread out this variations in this particular region. One could take this 
effect in account in an amplitude fit function using Monte Carlo studies. But the inclusion 
of phenomenological coherence factor could serve to mimic this effect and others possibles 
associated at the non-parametrized simple re-scattering contributions and the average on 
the isospin amplitudes. 

5 Closing Remarks 

Charmless three body B decays has an important appointment soon with the LHCb data 
[17]. In fact for the LHC run 2010-2012 it is expect more than decades of thousand events 
for the decays to three light mesons and thousand of events to three body decays 
channels involving a proton anti-proton [48] . This amount of events corresponds about 
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one order of magnitude bigger than the produced by the electron positron factories Belle 
and BaBar until today. With this statistics many new things could be understood for CP 
violation, that is the main proposal of LHCb experiment. Model independent approach 
was developed to search new sources of directly CP violation in Dalitz plot j49j . However 
if one want to go further and extract more information from Dalitz plot, like the CKM 
phase 7 [501 iSI] or the relative amount from the Penguins and Tree contributions or also 
to know the relative CP violation from each intermediary resonance state, or even to 
extract the the y parameter from 5° [52] one have to work with amplitude analysis. But 
as we could see in this paper, many studies have to be performed before to get trust-able 
parameters from this kind of fits. 

We could see here that the first results from Belle and BaBar presented some problems 
to describe the non-resonant contribution and the high mass scalar particles. The inclusion 
of an empirical exponential function increased significantly the confidence level of the fit. 
However we noted that this result imply in a strong destructive interference between the 
Si3 and S23 non-resonant components. This destructive interference take place just in the 
crossing region between the tt+tt" and K^tt~ resonances. We think one of the first studies 
to do with the new coming data is to look deeply at this effect. As it was pointed out in 
the previews section, the correct inclusion of the crossing region is the big challenger and 
seems that the charmless B three body decays analysis have to take care of this problem 
in order to have a good representation of the big amount of events expected for these final 
state B decays. 

The start point of the amplitude analysis for these studies could keep the general form 
of the isobaric model, but given more freedom to the fit function allowing more flexibility 
to understand better the data. Although Eq. 4 shown a smooth variation of the phase 
difference with data, this effect can be important for the fit involving a big amount of data 
as expected in LHCb and them could be add at the amplitude. In order to understand 
the importance of the coherence in hadronic interaction, it would be interesting includes 
a phenomenological real parameter in some interference terms of the fit function. The 
coherence parameters values, obtained from fit can give some light for understand better 
how works interferences in the crossing region among two resonances placed in different 
invariant mass, like tt+tt^ and K+vr" resonances in — ^ K^-k^ti^ decay. 

Acknowledgments: I would like to thanks Jussara Miranda, Cladio Lenz, Yakov 
Azimov and Alberto Reis for many suggestions and clarifying conversations about issues 
related with coherence in quantum mechanics. 
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